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1  Introduction 

In  our  ARO  work,  we  have  carried  out  an  extensive  research  program  for  the  study  of  robust  system  control 
using  methodologies  from  interpolation  theory,  dilation  theory,  and  functional  analysis.  We  have  also  become 
interested  in  image  processing  and  computer  vision,  and  their  application  to  visual  tracking  problems. 

First  of  all,  we  devoted  a  large  portion  of  our  research  effort  to  nonlinear  systems.  Thu  has  led  us  to 
derive  an  tieraltve  eommuiant  lifting  theorem  which  pves  an  explicit  design  procedure  for  nonlinear  systems 
and  captures  the  R"”-control  problem  in  the  nonlinear  framework  [38],  [39],  [40].  In  this  area,  we  have 
defined  a  notion  of  rsltonoltty  for  nonlinear  systems,  and  we  have  proven  that  the  iterative  commutant 
lifting  procedure  produces  rational  controllers  (in  this  nonlinear  sense)  if  we  start  from  rational  data.  The 
procedure  has  already  been  applied  to  certmn  systems  (connected  with  the  National  Aerospace  Plane) 
with  input  saturations  (in  collaboration  with  colleagues  at  the  Systems  Research  Center  of  Honeywell  in 
Minneapolis).  This  framework  has  also  led  to  new  directions  in  introducing  notions  of  causality  [41]  in 
commutant  lifting  theory,  and  moreover  has  led  to  the  possible  formulation  of  a  global  nonlinear  commutant 
lifting  theorem  as  a  sad^e-point  result.  We  are  particularly  interested  in  the  treatment  (and  understanding 
the  control  limitations)  of  "hard”  nonlineruities  such  as  dead-sone,  backlash,  and  saturation. 

During  the  research  period  supported  by  DAAL-03-91-G-0019,  we  also  developed  novel  interpolation 
methods  which  are  not  norm-based.  These  have  arisen  out  of  our  research  into  the  multivariable  gain  margin 
problem  as  well  as  the  more  general  itruetnred  eingnlar  value.  We  believe  that  this  new  type  of  interpolation 
theory  should  lead  to  several  interesting  directions  in  operator  theory  as  well. 

More  precisely,  our  work  on  the  multivariable  gain  margin  problem,  has  led  to  a  novel  interpolation  scheme 
which  we  call  spectral  Nevanliiuia-Pick  interpolation  and  the  more  general  structured  interpolation 
theory;  see  [14],  [13],  [15],  and  [17].  This  involves  matrix  interpolation  in  which  one  bounds  the  spectral 
radius,  and  not  the  norm  of  the  interpolants  as  in  ordinary  Nevanlinna-Pick  theory.  (Ordinary  Nevanlinna- 
Pick  is  precisely  the  type  of  mathematical  problem  that  arises  in  standard  synthesis.)  We  now  have  a 
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theoretical  procedure  for  solving  this  spectral  interpolation  problem  both  in  its  pure  matricial  and  tangential 
formulations. 

Our  solution  involves  a  generalisation  of  the  commutant  lifting  theorem,  and  a  new  object  of  interest 
in  linear  algebra  and  operator  theory  which  we  call  the  T-spectroi  radius.  We  are  presently  studying  both 
the  theoreticsd  and  practical  ramifications  of  our  solution,  and  we  would  like  to  develop  software  for  the 
computation  of  the  T*spectral  radius.  In  our  new  research,  we  have  been  concentrating  on  the  generalisation 
of  these  results  to  the  structured  singular  value  of  Doyle  and  Safonov  in  order  to  try  to  develop  an  analytic 
/s-synthesis  procedure  [15].  Closely  related  to  this  work  is  a  novel  lifting  technique  which  allows  us  to 
study  the  robust  stability  analysis  of  systems  under  various  kinds  of  structured  perturbations  which  we  will 
conuder  below.  (See  also  [9],  [10],  and  [05].)  This  Ufting  techmque  allows  one  to  interpret  the  upper  bound 
for  the  structured  singular  value  /i  defined  in  terms  of  cert^  scalings  as  a  structured  singular  value  on  a 
larger  space.  Moreover,  we  will  discuss  below  a  new  approach  for  the  approximation  of  Riemann  mappings 
(conformal  equivalences  from  simply  connected  sub-domuns  of  the  complex  plane  to  the  unit  disc)  for  the 
application  of  solving  the  g^-phase  margin  problem  [48]. 

We  have  continued  our  research  on  the  utilisation  operator  theoretic  methods  in  IT'^-optimisation  theory 
using  skew  ToepUix  operuiort  (see  [12]  for  the  precise  definition),  which  seem  ideally  suited  for  studying 
design  problems,  especially  for  d^tributed  (i.e.,  infinite  dimensional)  systems.  These  methods  are  quite 
natural  in  the  control  context  since  they  allow  one  to  do  design  just  using  the  input/output  operators. 
This  has  led  to  an  explicit  solution  of  the  four  block  problem  for  large  class  of  multivariable  distributed 
systems.  (This  is  the  most  general  R°°-optimisation  problem.)  The  procedure  seems  to  be  numerically 
robust  as  evidenced  by  its  implementation  at  the  Systems  Research  Center  of  Honeywell,  and  its  application 
to  a  number  of  distributed  plant  models  including  a  flexible  beam  [52],  and  an  unstable  delay  system 
associated  with  the  flight  control  of  the  X>29  [26].  A  nice  feature  of  this  approach  is  that  the  complexity 
of  the  computations  only  depends  on  the  the  weighting  matrices  (modelling  the  disturbances)  and  not  on 
the  plant  (which  may  be  distributed).  Since  the  weighting  matrices  are  typically  taken  to  be  rational,  this 
approach  seems  very  efficient  even  in  the  finite  dimensional  case  for  plants  with  large  state  spaces.  We  have 
also  been  developing  methods  for  the  utilisation  of  an  operator  due  to  N.  Young  in  order  to  simplify  some 
of  our  computations  in  the  ttandard  problem  for  multivariable  systems. 

Next  using  the  one-step  dilation  procedure,  we  have  given  a  way  of  parametrising  the  suboptimal  con¬ 
trollers  for  such  generalised  interpolation  problems  [34].  This  has  led  to  a  scheme  for  designing  suboptimal 
finite  dimensional  controllers  for  cflstributed  systems  in  the  one  block  setting  [57].  We  would  like  to  con¬ 
tinue  this  direction  for  the  general  standard  problem  as  well  as  apply  the  techniques  to  a  number  of  design 
problems.  Some  steps  in  this  direction  have  already  been  taken  in  [26]  and  [52]. 

Closely  related  to  the  above  work  has  been  the  use  of  techniques  to  study  sampled-data  control. 
See  [7]  and  the  references  therein.  Here  one  can  use  a  certmn  operator-theoretic  lifting  method  to  describe 
a  complete  solution  to  the  analysis  problem  of  verifying  that  a  given  controller  constrains  the  X’-induced 
norm  of  the  sampled-data  system  to  be  less  than  some  pre-specified  level.  (In  fact,  the  lifting  method  is 
applicable  to  all  norm  based  optimisation  problems.  In  discrete-time  such  a  technique  was  employed  in  [46] 
to  study  a  number  of  issues  in  robust  control.) 

We  will  also  discuss  below  our  research  in  combined  suboptimal  controllers  ([27]  and  [29])  as 

well  as  sketch  an  approach  for  combining  other  norms  based  on  techniques  from  interpolation  on  Banach 
spaces,  and  the  connection  of  this  methodology  to  a  new  approach  for  model  reduction  which  we  intend  to 
explore  in  our  new  ARO  contract. 

In  a  different  direction,  we  have  been  using  curve  evolution  theory  (and  especially  invariant  nonlinear 
diffusion  equations),  for  the  development  of  new  algorithms  in  image  processing  and  computer  vision.  Curve 
evolution  theory  has  recently  become  a  major  topic  of  research,  and  indeed  has  appeared  in  such  diverse 
fields  as  differential  geometry,  parabolic  equations  theory,  numerical  analysis,  computer  vision,  the  viscosity 
solutions  of  Hamilton-Jacobi  equations,  and  image  processing.  See  [49,  50,  51,  60,  61,  62,  63,  64]  and  the 
references  thernn. 

In  particular,  evolution  equatioiu  which  are  geometric  non-linear  versions  of  the  classical  heat  equation 
have  received  much  attention,  since  these  equations  have  both  theoretical  and  practical  importance.  These 
ideas  have  been  used  also  in  shape-from-shading,  image  smoothing  and  enhancement,  motion  planning,  shape 
segmentation,  optical  flow,  and  a  continuous  implementation  of  mathematical  morphology. 
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A  complete  list  of  publications  which  acknowledges  support  of  DAAL-03-91-G-0019  is  included  in  Sec¬ 
tion  7.  We  now  sketch  the  work  done  by  A.  Tannenbaum  and  collaborators  on  this  contract. 


2  Nonlinear  Robust  Control 

We  have  been  pursuing  several  directions  in  order  to  derive  nonlinear  generalisations  of  the  (linear)  S°° 
theory  in  the  wdghted  sensitivity  minimisation  (one  block),  mixed  sensitivity  (two  block),  and  even  standard 
problem  frameworks.  One  such  direction  is  based  on  an  iterative  commutant  lifting  theorem  [38],  [39], 
[40],  [32]  which  pves  an  explicit  deugn  procedure  for  nonlinear  systems  and  captures  the  JJ°°-control  problem 
for  a  la^e  class  of  nonlinear  plants. 

2.1  Causal  Analytic  Mappings 

We  conuder  analytic  mappings  on  Hilbert  space.  For  the  precise  definitions  see  [38],  [39],  [40],  [30].  For 
simplicity,  we  will  only  consider  SISO  operators  in  what  follows.  The  multivariable  case  works  the  same  way 
[30].  We  will  find  it  convenient  to  employ  the  Fourier  representation  of  elements  of 

So  we  consider  an  analytic  map  ^  with  Q  —"K  (the  standard  Hardy  space  on  the  disc  D).  Note 

that 

®  ®  if*  =  (if’)*"  S 

where  we  map  (sin  the  i-th  place)  to  Zi,  t  =  1,  •  •  * ,  n.  In  the  usual  way,  we  say  that  ^ 

Mhifl-invariani  (or  time-invariant)  if 

*n5«"=S^„  Vn>l, 

where  S  :  denotes  the  canonical  unilateral  right  shift.  (Equivalently,  this  means  that  =  ^oS 

on  some  open  ball  about  the  oripn  in  which  ^  is  defined.) 

Now  set 

*“  ®  ®  P(j)  (n  times),  i  >  1,  »  >  1, 

where 

Then  we  say  that  ^  is  causal  if 

^  ^  I- 

For  ^  —*  linear  and  time-invariant,  it  is  easy  to  see  that  ^  is  causal.  In  the  nonlinear  setting 

however,  time-invariance  may  not  imply  causality  [40].  It  is  for  this  reason,  that  a  causality  constraint  must 
be  explicitly  included  for  nonlinear  design. 

2.2  Causal  Optimization  Problem 

Because  of  the  difficulties  involving  causality  when  one  applies  the  classical  commutant  lifting  theorem  in 
the  nonlinear  framework,  we  will  need  to  formulate  a  new  linear  causal  optimization  problem.  Then 
we  will  indicate  how  to  reduce  the  nonlinear  generalisation  of  the  H'"  sensitivity  minimisation  problem  to 
a  series  of  such  problems. 

We  let  5(n)  denote  the  unilateral  shift  on  given  by  multiplication  by  (zi  •  -  -Sn).  Since 

will  be  fixed  in  the  discusnon  we  will  let  5  :=  5(„).  In  what  follows,  U  will  denote  the  unilateral  shift  on 
given  by  multiplication  by  z,  and  0  €  JT**  wffi  be  an  inner  function.  Finally  W  :  will 

denote  a  causal,  time-invariant  bounded  linear  operator. 

We  can  now  state  the  eanzal  -optimization  problem  (COP):  Find 

tr  :=  inf{||W  -  0Q|| :  Q  :  Jr*(P”)  — ♦  H*,  Q  causal,  time-invariant}.  (1) 

Moreover,  we  want  to  compute  an  optimal,  causal,  time-invariant  Qgpt  such  that 

o  =  |iw-eQ,,||.  (2) 
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Figure  1:  Sensitivity  Minimisation  Configuration 


If  we  drop  the  causality  constraint,  the  solution  to  problem  (1)  is  provided  by  the  classical  commutant 
lifting  theorem.  With  the  causality  construnt,  the  solution  to  (COP)  is  abstractly  provided  by  a  causal 
commutant  lifting  theorem  [41],  [33],  [31]. 

Our  constructions  are  baaed  on  a  rtdueiion  tkeortm  [30],  which  allows  us  to  reduce  the  causal  J7°°- 
optimisation  problem  to  one  of  classical  interpolation.  Then  in  analogy  with  our  previous  work,  we  can  base 
an  causal  iterative  design  procedure  on  this  chiun  of  ideas. 

Here  is  how  such  a  procedure  would  look.  Let  us  call  an  analytic  input/output  operator  ^ 
admuiibU  if  is  is  causal,  time-invariant,  and  ^(0)  =  0.  Denote  the  set  of  admissible  operators  by  Ca.  In 
what  follows  below,  we  assume  P,  W  €  Ca,  and  that  W  admits  an  admissible  inverse. 

Referring  to  Figure  1,  we  consider  the  problem  of  finding 

fit  :=  inf  sup  ||[(J  +  P  o  C)“^  o  W]w||,  (3) 

®  IWIS* 

where  we  take  the  infimum  over  all  stabilising  controllers.  (In  what  follows,  we  let  ||  ||  denote  the  2-norm 
II  IIj  on  as  well  as  the  assodated  operator  norm.  The  context  will  make  the  meaning  clear.)  Thus  we  are 
looking  at  a  worst  case  disturbance  attenuation  problem  where  the  energy  of  the  signals  v  is  required  to  be 
bounded  by  some  pre-spedfied  level  6.  (In  the  linear  case  of  course  since  everything  scales,  we  can  dways 
without  loss  of  generality  take  f  =  1.  For  nonlinear  systems,  we  must  specify  the  energy  boimd  a  priori) 
Then  one  sees  that  (3)  is  equivalent  to  the  problem  of  finding  the  problem  of  finding 

pi=  W  sup  ||(W-Po9)t>||.  (4) 

The  iterated  causal  commutant  lifting  procedure  gives  an  approach  for  approximating  a  solution  to  such 
a  problem.  Briefly,  the  idea  is  that  we  write 

W  =  Wi  +  W2  +  -'-, 

P  =  Pi  +  Pa  +  •  •  • , 

q  =  qi+ga-i - , 

where  Wj,Pj,qj  are  homogeneous  polynomials  of  degree  j.  Notice  that 

fit  =  6  M  \\Wi  -  Pi, ill  +  0(«>),  (5) 
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where  the  latter  norm  U  the  operator  norm  (i.e.,  norm).  From  the  classical  commutant  lifting  theorem 
we  can  find  an  optimal  (linear,  causal,  time-invariant)  qi,api  €  H‘“  such  that 

W=tf||Wi-PHi..pt||  +  0(tf»).  (6) 

Now  the  iterative  procedure  gives  a  way  of  finding  higher  order  corrections  to  this  linearisation.  Let  us 
illustrate  this  now  with  the  second  order  correction.  Indeed,  having  fixed  the  linear  part  91, opt  of  q  in  (4), 
we  note  that 

W(i,)  -  P{q{v))  -  {Wi  -  Piqi.„t){v)  = 

Wa(«)  -  ft(«i.opt(i»))  -  Piqiiv)  +  higher  order  terms. 

Regarding  W2,  l^,9a  »  linear  operators  on  SL  C)  as  above,  we  see  that 

sup  \\{W  -Po  q){v)  -  {Wt  -  Pxqi.^i)v\\  <  6^\\W2  -  Px,al|  0(4’), 
where  the  “weight"  W2  is  given  by 

W2  :=  W2  —  Pa(9i,op(  ®  9x,opf)- 

The  point  of  the  iterative  causal  commutant  lifting  procedure  is  to  allow  us  to  construct  an  optimal  admissible 
qijopu  eixd  so  on. 

In  short,  instead  of  simply  designing  a  linear  compensator  for  a  linearisation  of  the  given  nonlinear  system, 
this  methodology  allows  one  to  explicitly  take  into  account  the  higher  order  terms  of  the  nonlinear  plant,  and 
therefore  increase  the  ball  of  operation  for  the  nonlinear  controller.  Moreover,  if  the  linear  part  of  the  plant 
is  rational,  our  iterative  procedure  may  be  reduced  to  a  series  of  finite  dimensional  matrix  computations. 
(See  [30,  38]  for  discussions  of  rationality  in  the  nonlinear  framework.) 

Finally,  in  a  number  of  interesting  cases,  the  above  methodology  can  be  extended  to  a  nonlinear  extension 
of  the  mixed  sensitivity  problem;  see  [25].  In  fact,  one  can  show  that  for  linear  weighting  filters  and 
adnussible  plant,  the  nonlinear  mixed  sensitivity  problem  may  be  reduced  to  a  standard  linear  two  block 
problem,  followed  by  a  one  block  nonlinear  design,  which  may  be  solved  using  the  above  iterative  methods. 

Recently  the  above  procedure  was  extended  to  the  general  standard  prohlem  for  nonlinear  systemr,  see 
[32].  We  should  note  that  we  have  concentrated  on  time-invariant  systems  in  our  above  description.  In  fact, 
there  is  a  much  more  general  causal  commutant  lifting  methodology  [41,  33,  31]  that  can  be  extended  in 
principle  to  time-varying  systems.  This  is  a  very  important  problem  area.  Moreover,  there  are  interesting 
problems  concerning  the  equality  of  certain  invariants  which  arise  the  the  general  causal  commutant  lifting 
framework  [41]  which  we  intend  to  consider  in  our  upcoming  ARO  contract. 

2.3  Saddle-Point  Approach  to  Nonlinear  Optimization 

We  have  discussed  above  an  iterative  commutant  lifting  approach  to  nonlinear  system  design.  The  iterative 
communtant  lifting  technique  is  basically  a  local  analytic  method  for  nonlinear  system  synthesis.  We  have 
also  been  exploring  a  very  different  global  approach  based  on  an  interpretation  of  the  classical  commutant 
lifting  theorem  as  a  saddle-point  result  [42].  This  motivates  us  to  formulate  a  nonlinear  commutant  lifting 
result  in  such  a  saddle-point,  game-theoretic  framework.  A  related  approach  to  nonlinear  design  has  already 
been  independently  employed  in  the  novel  and  important  papers  of  Ball-Helton  [4],  [5]. 

In  our  research,  instead  of  considering  general  nonlinear  systems  we  limit  ourselves  to  the  concrete  (but 
certainly  interesting  case)  of  linear  systems  with  input  saturations.  Such  systems  occur,  of  course,  all  the 
time  in  “nature.”  We  should  add  that  a  similar  approach  should  be  valid  for  many  of  the  hard,  memoryless, 
noninvertible  nonlinearities  which  appear  in  control. 

Mathematically,  the  case  of  the  saturation  is  very  interesting  since  in  a  certun  sense  it  is  a  nonlinear 
analogue  of  an  inner  (non-minimum  phase)  element,  whose  “spectrum”  seems  to  be  spread  throughout 
the  unit  circle.  Thus  the  problem  of  sensitivity  minimisation  for  such  elements  (that  is,  a  form  of  weighted 
inversion)  is  particularly  dallenging,  and  will  be  a  key  topic  to  be  studied  in  our  upcoming  research  program. 
What  needs  to  be  developed  for  this  global  approach  to  nonlinear  is  a  commutant  lifting  type  result 
valid  on  convex  spaces.  We  are  doing  this  in  our  new  ARO  contract. 
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3  Structured  Singular  Values  and  ;x-Synthesis 

We  would  like  now  to  diacuss  our  work  on  the  structured  singular  value  and  /^synthesis  using  a  novel 
approach  which  we  developed  called  ttnclurti  interpolation  theory.  The  structured  singular  value  was 
introduced  into  robust  control  by  John  Doyle  and  Michael  Safonov  ([22],  [59])  to  handle  problems  involving 
structured  perturbations  which  includes  both  and  the  multivariable  gun  margin  as  special  cases. 

Our  starting  point  is  the  fact  that  the  problem  of  internal  stabilisation  of  a  given  LTI  multivariable  plant 
r*n  be  reduced  to  one  of  interpolation.  Thus  the  natural  measure  of  robust  performance  in  this  framework 
is  given  by  the  minimisation  of  the  structured  singular  values  over  all  possible  interpolants.  Hence,  one 
needs  a  structured  Nevanlinna-Pick  type  result  which  will  generalise  both  the  classical  Nevanlinna^Pick 
theorem  (relative  to  the  H°°-norm),  as  well  as  the  spectral  Nevanlinna-Pick  theorem  (relative  to  the  spectral 
radius).  Such  a  structured  Nevanlinna-Pick  theorem  can  be  deduced  as  a  consequence  of  a  general  structured 
commutant  lifting  theorem  as  we  will  indicate  below.  This  approach  has  been  developed  in  our  ARO 
sponsored  work  [13,  14,  11,  17,  8]. 

We  will  discuss  some  of  these  issues  now,  and  explun  how  these  results  can  serve  as  a  basis  of  performing 
the  /i-synthesis  procedure  in  robust  feedback  control  in  a  rigorous,  analytical  manner. 

3.1  Structured  Singular  Values  and  Dilations 

We  would  like  to  formally  introduce  the  structured  singular  value  now,  and  give  some  of  its  basic  properties. 
We  base  this  discussion  on  [IS].  Instead  of  working  over  diagonal  sets  of  matrices  as  in  [22],  we  can  more 
generally  work  over  an  arbitrary  finite  dimensional  C*-algebra. 

Let  £  be  a  complex  finite-dimensional  Hilbert  space,  and  A  C  jC(f )  (the  space  of  bounded  linear  operators 
on  €),  a  C*-algebra.  For  A  €  £(£),  A^0,ve  define  the  structured  singular  value 

Ma(A)  :=  pnf{||A'|| :  JT  €  A.  -1  €  o'(AJr)}]-»  . 

Moreover,  we  set 

M^)  :=inf{|lJrAA-M|:Jr€A'}, 

where  A'  is  the  commutator  of  A.  Note  that  for  A  =  jC(S),  pa(A)  =  ||A||,  while  for  A  =  C/f ,  A>a(A)  =  ||A||(p 
(the  spectral  radius  of  A). 

For  certain  diagonal  algebras  of  matrices,  it  is  argued  in  [22]  and  [59]  that  the  structured  singular  value 
Pa  is  the  natural  object  of  study  in  robust  control.  Unfortunately  the  structured  singular  value  is  difficult 
to  compute,  so  in  practice  it  is  pa  which  is  actually  used  for  control  problems.  It  is  therefore  of  interest  to 
know  when  these  two  objects  are  equal. 

In  [22],  Doyle  has  shown  that  in  fact  pa  =  Pa  when  the  relevant  diagonal  algebra  has  three  or  fewer 
blocks.  In  [15],  we  give  a  very  different  proof  of  this  fact  based  on  the  following  result  which  we  believe  has 
independent  interest.  Define  the  operator  Ma  €  £(^)  by  Ma  ~  AX.  Notice  that  C(€)  may  be  given  a 
Hilbert  space  structure  with  respect  to 

(Ti,T3)  :=rr(T;Ti), 

where  Tr  denotes  the  trace.  Define 

Pa(A):=  PA(3f^) 

where 

A  :=  {Mx  :  X  €  A'}'. 

We  now  have  (see  [15]  for  the  proof): 

Theorem  1 

PA  (A)  :=  Pa  (A). 
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Note  that  Theorem  1  implies  that  Ji^  can  be  regarded  as  a  structured  singular  value  on  a  big' 
ger  space.  From  the  theorem  and  property  (ii)  above,  we  can  immediately  infer  the  key  fact  that  is 
continuous.  This  result  is  also  strongly  connected  to  some  recent  work  on  robust  stability  with  respect 
to  time-varying  perturbations  [65].  In  fact,  we  are  investigating  the  extension  of  this  lifting  technique  to 
operators  on  infinite  dimensional  Hilbert  spaces  with  applications  to  the  robust  stability  analysis  of  systems 
under  various  kinds  of  structured  perturbations  [9].  In  particular,  we  have  considered  the  following  two  cases 
[9,  10]: 

(i)  The  algebra  of  constant  diagonal  scales  A',  with  the  operator  A  taken  to  be  analytic  Toeplits.  (This 
is  the  case  analysed  in  [65].)  Note  that  by  “analytic  Toeplits  operator”  A,  we  mean  that  A  is  given 
by  an  n  X  n  matrix  with  entries  If*”  functions,  which  we  regard  as  acting  as  a  multiplication  operator 
on  J7’(C”).  Thus  A  defines  a  stable  LTI  system.  We  also  want  to  consider  general  time-varying  linear 
input-output  operators  with  this  constant  (Uagonal  algebra  of  scalings  as  well. 

(ii)  Ag^  we  take  A  to  be  analytic  Toeplits,  but  now  we  want  to  consider  A'  to  be  the  algebra  of  diagonal 
analytic  Toeplits  operators.  This  is  the  type  of  problem  considered  in  p-synthesis. 

We  should  note  that  this  work  has  also  led  to  a  new  relative  ToepliU-Hautdorff  theorem. 

We  will  now  outline  a  structured  analogue  of  the  commutant  lifting  theorem  [15].  This  will  be  applied 
to  the  structured  version  of  classical  matricial  Nevanlinna-Pick  interpolation  below. 

Set  T  :=  5(m)  ®  Je  where  m  is  a  finite  Blaschke  product,  5(m)  is  the  compressed  shift,  and  £  is  a  finite 
dimensional  complex  Hilbert  space.  Fix  A  C  £(£),  a  C*-algebra.  Define 

Iffi  0  A  :=  {//r>  ®  :  Jf  €  A},  Ijfi  ®  A‘  :=  {/^a  ®  :  X  6  A*}. 

Notice  that  H  :=  0  mH^{£)  reduces  both  ®  A  and  Igt  ®  A'.  Now  define  for  A  €  {T}'  (the 

commutant  of  T), 

^(A)  :=  inf{l|jrAX-»|| :  X  invertible,  X  6  {T}',  X  6  (/»»  0  Al£)'}. 

Let  H  be  the  isometric  dilation  of  T  on  £r^(£)  (so  that  IT  is  defined  by  multiplication  by  z).  ForH  €  {C''}', 
define 

p%{B)  :=  inf{|iyHy-^||  :  Y  invertible,  Y  G  Y  G  (/»*  0  A)'}. 


We  can  now  state  the  following  result  from  [15]: 


Theorem  2  (Structured  Commutant  Lifting  Theorem)  Notation  as  above.  Then  for  A  G  {T}', 

Pr(A)  =  inf{p^(H)  :  B  it  a  commuting  dilation  of  A}. 


Actually,  in  the  finite  dimensional  case  of  interest  to  us,  we  can  show  that 

/^(A)  =  irf  sup  iiA(B(z)) 

where  B  G  {U}'  is  a  commuting  rational  dilation  of  A.  This  leads  to  structured  versions  of  the  classical 
Nevanliima-Pid  theory  both  in  its  matricial  and  tangential  forms  [15, 11,  8].  The  tangential  form  is  precisely 
what  is  needed  in  making  the  ^-synthesis  procedure  rigorous. 
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3.2  Conformal  Mappings 

We  would  like  to  briefly  discun  some  work  that  we  also  did  on  the  gain-phase  margin  problem  that  involves 
some  novel  approximation  methods  of  Riemann  mappings  (conformal  equivalences  from  a  simply  connected 
proper  subset  of  C  to  the  unit  disc)  due  to  Marshall-Morrow  (see  [20]  and  the  references  therein). 

The  successful  solution  of  the  gain  and  phase  margin  problems  (see  [67],  [48],  [23])  was  based  on  the  fact 
that  one  could  explicitly  construct  a  conformal  equiv^ence  from  a  simply  connected  subdomain  of  C  (which 
was  associated  to  the  uncertunty  model)  to  the  disc.  Nevanlinna-Pick  interpolation  then  took  care  of  the 
rest.  In  [48],  we  considered  a  combination  of  both  the  classical  gain  and  phase  margins  in  order  to  arrive  at  a 
better  measure  of  robust  stability  which  we  called  the  yata-pkose  margin.  Unfortunately,  in  this  case,  there  is 
no  explicitly  computable  conformal  equivalence  between  the  associated  domain  of  uncertainty  and  the  disc. 
The  algorithm  of  Marshall-Morrow  however  gives  a  fast  way  of  approximating  the  given  equivalence,  and 
has  led  to  a  reliable  approcdmate  solution  of  gain-phase  margin  in  [20].  We  are  now  applying  this  algorithm 
to  the  type  of  general  robust  synthesis  question  as  considered  by  Sideris  and  Safonov  in  [66]. 

4  H°°  Optimization  of  Distributed  Systems 

In  [53],  [54]  we  have  considered  the  mixed  sensitivity  j?‘"-optimisation  problem  for  distributed  plants  with 
a  finite  number  of  unstable  poles. 

As  in  the  standard  operator  approach  to  such  problems,  the  computation  of  the  optimal  performance 
and  corresponding  optimal  controller  is  reduced  to  a  finite  dimensional  matrix  problem.  In  the  stable  case, 
the  sise  of  the  matrix  only  depends  on  the  McMillan  degree  of  the  weighting  filters.  In  the  case  of  unstable 
plants,  the  sise  of  the  corresponding  matrix  depends  on  the  number  of  right  half  plane  poles  of  the  plant 
as  well.  The  dimension  of  this  matrix  can  be  computed  a  priori.  The  key  mathematical  fact  that  we  use  is 
that  the  skew  Toeplits  operators  which  we  obtain  in  the  unstable  case  are  finite  rank  perturbations  of  the 
classical  skew  Toeplits  operators  obtained  from  compressions  of  rational  functions. 

In  our  ARO  research,  we  have  used  these  techniques  to  synthesise  controllers  for  several  types  of  flexible 
structures  [52],  and  an  unstable  delay  system  derived  from  the  flight  control  of  the  X-29  [26]. 

4.1  Two  Block  Problems  With  Unstable  Plant 

In  this  section,  we  will  show  that  several  two  block  Jf”-minimisation  problems  reduce  to  the  computation 
of  the  norm  of  a  certain  skew  Toeplits  operator,  and  indicate  how  this  norm  may  be  computed.  We  begin 
with  some  notation.  The  Hardy  spaces  and  are  defined  on  the  unit  disc  as  above.  We  denote 

RH°°  :=  {rational  functions  in  H°°}. 

We  consider  the  feedback  configuration  of  Figure  2  with 


and  Gn  €  Gi  e  RH^. 

We  assume  that  (i)  G„  =  m„Gn«,  where  m„  €  is  inner  (arbitrary)  and  G^q  €  is  outer,  and  (ii) 
Gn,  Gi  have  no  common  seros  in  the  closed  unit  disc.  We  also  write  Gi  =  miGie  where  mj  €  is 

inner  and  G*,  €  AH*  is  outer.  Under  these  assumptions  there  exist  X  €  AH"  and  y  €  H"  such  that 

HG„  +  yGd=l.  (7) 

The  set  of  all  controllers  which  stabilise  the  plant  can  now  be  written  in  the  form 
X  +  QGj 
Y-QGn 

for  some  Q  €  H".  Now  let  S  :=  (H-  PC)~^  and  note  that 

S=l-XGn-QGnGi.  (8) 
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Figure  2:  Standard  Feedback  Configuration 


In  [53],  [54],  we  show  that  the  computation  of 

=  ...wSUc  I  [  Wj(S  -  1)  ]  1 

where  IFi,  IFa  €  are  given  weighting  functions  with  6  RH°°  may  be  reduced  to  computing 

the  norm  of  the  operator 

A  .=  (Wo(S)  -  Wo{S)m{S)) 

Go(S)  J  ’ 

where  S:  —*  denotes  the  unilateral  shift,  H{my)  :=  H^Qm^  and  Pj7(m,)  the  orthogonal  projection 

onto  fir(ni«),  for  m,m«  inner  functions  associated  to  the  plant  and  weighting  filters,  and  where  Wo,H'b«Go 
are  rational  fiT"’  functioiu  computed  from  the  plant  and  weighting  filters.  This  reduction  is  true  for  plants 
with  arbitrary  outer  parts. 

In  [53],  [54],  we  develop  an  approach  to  computing  the  singular  values  and  vectors  of  operators  of  the 
form  (9).  We  remark  here  that  it  is  easy  to  compute  the  essential  norm  (see  [53],  [54])  of  the  operator  A, 
which  will  be  denoted  by  ||A||,.  Then  in  [54],  the  following  result  is  proven; 


Theorem  3  Let  n  denote  the  maximum  of  the  McMillan  degrees  of  the  weighting  filters  Wi  and  Wj,  and 
let  t  denote  the  number  of  unstable  poles  of  the  plant  P.  Then  the  singular  vectors  and  valves  of  A  which 
are  >  ||A||«  may  be  derived  from  an  esplicHly  computable  system  o/3n  +  2/  linear  equations  (the  “singular 
system"). 


In  [54],  the  singular  system  of  equations  is  explicitly  written  down.  Again  the  number  of  equations 
only  depends  on  the  McMillan  degrees  of  the  weighting  filters,  and  the  number  of  right  half  plane  poles  of 
the  plant.  The  computation  of  the  maximal  singular  value  and  the  associated  singular  vectors  of  A,  then 
allows  us  to  find  the  optimal  performance  p  of  our  original  control  problem  and  the  corresponding  optim2d 
compensator. 

4.2  Mixed  Norm  Suboptimal  Controllers 

One  of  the  main  themes  in  modern  control  theory  is  the  utilisation  of  norm  based  criteria  to  measure  the 
optimal  performance  of  a  pven  control  system.  (See  [23]  for  detuled  discussions  on  this  topic.)  Two  of  the 
most  important  norms  employed  in  modem  control  analysis  and  design  are  the  ff^-norm  and  the  J7°°-norm. 
Indeed,  the  ff^-aorm  is  the  basis  of  classical  quadratic  optimal  control,  while  the  IT*” -norm  appears  in 
modem  robust  synthesis  (and  implicitly  in  the  more  classical  loop  shaping  methods  [23]). 

We  have  been  investigating  combining  the  advantages  of  both  and  control  in  [29].  (See  also  [27] 
for  extensive  lists  of  references  on  various  approach^  to  mixed  norm  control.)  Our  starting  point  is  a  nice 
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result  of  Kftftal-Larson-Weiss  [45]  which  guarantees  the  existence  of  an  interpolant  which  simultaneoiuly 
satisfies  an  and  suboptimality  criterion.  Using  the  theory  described  in  [34]  on  suboptimal 
interpolants,  we  can  then  give  an  explicit  way  of  computing  the  combined  suboptimal  interpolant 

in  a  given  f**  problem. 

More  precisely,  for  c  €  L’”,  set  dao(c)  :=  di8tance(c,  H°°),  and  d2(c)  :=  di8tance(c,  H^).  (Note  that  all  of 
our  Hardy  spaces  are  defined  on  the  unit  disc  here.)  Then  for  any  6  >  1,  it  is  proven  in  [45]  that  there  exists 
^  €  i?*”  su^  that 

|lc-^||co<Woo(c),  ||c-^|j2<-^^.  (10) 

In  case,  m^w  €  with  m  inner,  and  vt  =  p/q  rational,  using  skew  Toeplits  theory  [29]  one  may 
explicitly  construct  i-e-i  if  e  :=  mw,  then  we  can  find  ^  satisfying 

||mto  -  ^lloo  <  Woo(c),  ||mto  -  ^lls  <  (11) 

Indeed,  let  p  :=  Sdoo(c),  and  let  Ap  :=  P^q(Tyq(T)  —  p(T)*p(T)  be  a  skew  Toeplits  operator.  Let  q/  denote 
the  function  in  defined  by 

qt  :=  6®doo(c)*i~^«(T)*(l  -  Tnm(0)).  (12) 


Then  the  function 


w  —  Tn4  =  Bt 


p(T)qi 

mm(0)  +  q{T)qt 


(13) 


satisfies  (11). 

The  skew  Toeplits  approtich  is  the  only  procedure  of  which  we  are  aware  to  compute  a  ^  satisfying  (11) 
when  m  is  irrational.  Further,  let  n  :=  max{degree  p,  degree  q).  Then,  if  m  is  Blaschke  product  of  order  k 
where  n  is  small  and  k  is  large,  the  methodology  just  outlined  provides  a  numerically  efficient  method  for 
computing  the  function  ^  satisfying  the  Kaftal-Larson-Weiss  constraints  (11).  However,  if  this  is  not  the 
case  and  m  is  rational,  then  the  state  space  methods  in  [27]  to  compute  a  ^  satisfying  (11)  may  be  more 
efficient. 

Since  as  it  turns  out  that  the  combined  interpolant  is  a  central  solution,  these  interpolants  are 

well-known  to  be  numerically  robust,  and  leads  to  an  interesting  class  of  combined  controllers  for 

feedback  systems.  We  are  now  actively  investigating  the  properties  of  these  controllers,  as  well  as  extending 
these  results  to  find  a  computational  procedure  for  the  fidl  standard  problem  setting  for  distributed  MIMO 
systems. 

Further,  we  have  found  a  new  method  for  constructing  the  central  solutions  as  related  to  the  general 
standard  problem  in  the  multivariable  case.  We  are  now  in  the  process  of  exploiting  this  technique  for  the 
explicit  parametrisation  of  the  suboptimal  controllers  for  multivariable  distributed  systems.  This  would  give 
a  powerful  alternative  method  to  the  one-step  extension  techniques  we  have  been  previously  using. 

Finally,  we  are  interested  in  combining  other  types  of  norms  in  control,  for  example,  and  H°°.  A 
possible  method  for  doing  this  can  be  based  on  some  new  results  on  interpolation  in  Banach  spaces  due 
to  Pisier  [58].  In  fact,  this  methodology  can  also  have  some  nice  applications  to  model  reduction  (relative  to 
**best”  Hankel  norm  approximations),  which  is  certunly  a  direction  which  we  intend  to  explore. 


5  Image  Processing  and  Computer  Vision 

We  have  been  doing  extensive  research  into  image  processing  and  computer  vision,  and  the  use  of  these 
methodologies  in  visual  tracking.  We  are  very  interested  in  understanding  how  to  use  visual  information  in  a 
feedback  loop.  In  order  to  give  the  reader  a  flavor  of  this  work,  we  will  concentrate  on  a  new  multiresolution 
representation  theory  that  we  developed  for  planar  shapes  [60,  61].  In  this  work,  we  combine  the  three 
important  theories  of  scale-spaces,  affine  invariant  descriptors,  and  curve  evolution,  in  order  to  define  a 
new  aflSne  invariant  scale-space  for  plane  curves.  This  scale-space  is  obtmned  from  the  solution  of  a  novel 
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curve  flow,  which  admits  affine  invariant  solutions.  Our  other  work  in  vision  covers  a  new  theory  of  shape 
[49,  SO,  51],  as  well  as  new  methods  for  affine  invariant  image  processing  and  analysis  [62,  63,  64]. 

Multi-scale  descriptions  of  signals  have  been  studied  for  several  years  already.  A  possible  formalism  for 
this  topic  comes  from  the  idea  of  multi-scale  Altering  which  was  introduced  by  Witkin  [68],  and  developed  in 
several  diflerent  frameworks  by  a  number  of  researchers  in  the  past  decade.  The  idea  of  scale-space  Altering 
is  very  simple  and  can  be  formulated  as  follows:  Given  an  initial  signal  :  R-**  -*  R-'”i  tbe  scale-space 

is  obtained  by  Altering  it  with  a  kernel  K{X,i)  :  —*  R*",  where  t  €  R***  represents  the  scale.  In  other 

words,  the  scale-space  u  given  by  deAned  as 

♦(.?,l):=n^(^_,,[«o(^)],  (14) 

where  njc(.,t)[']  represents  the  action  of  the  Alter  AC(*,i).  Larger  values  of  t  correspond  to  images  at  coarser 
resolutions. 

It  is  important  to  note  that  not  every  kernel  can  be  used  in  deAning  a  scale-space.  Indeed,  several 
conditions  need  to  be  imposed  on  the  signal  9{X,t)  (and  therefore  on  the  Altering  operation  (14)).  One  of 
the  most  important  is  that  of  eamality,  which  states  that  no  "information”  is  created  when  moving  from 
Ane  to  coarse  scales. 

A  famous  example  of  a  kernel  which  satisAes  the  required  conditions  is  the  Gaussian  kernel.  In  this  case, 
the  scale-space  is  linear,  and  the  Alter  in  (14)  is  just  deflned  by  convolution.  The  Gaussian  kernel  is  one  of 
the  most  studied  in  the  theory  of  scale-spaces.  It  has  some  very  interesting  properties,  one  of  them  being 
that  the  rignal  '9  obtmned  from  it,  is  the  solution  of  the  heat  equation  (with  'fo  ns  initial  condition)  given 
by 


(For  more  details  about  the  Gaussian  scale-space,  see  the  aforementioned  references.) 

One  of  the  key  facts  that  can  be  gleaned  from  the  Gaussian  example,  is  that  the  scale-space  can  be 
obtained  as  the  solution  of  a  partial  differential  equation  called  an  evolution  eguafton.  This  idea  was  developed 
in  different  works  for  evolution  equations  different  from  the  classical  heat  equation.  In  what  follows,  we  will 
describe  scale-spaces  for  planar  curves  which  are  obtained  as  solutions  of  nonlinear  evolution  equations. 

The  second  fundamental  concept  which  we  would  like  to  emphasise  is  that  of  tnvamnf  deaeriptor.  An 
invariant  descriptor  is  a  property  of  an  object,  which  does  not  change  when  the  object  undergoes  certain 
transformations.  More  precisely,  a  quantity  Q  is  called  an  invariant  of  a  Lie  group  Q  if  whenever  Q  transforms 
into  Q  by  any  transformation  g  €  Q,  yte  obtun  Q  =  QQ,  where  6  is  a  function  of  g  alone.  If  6  =  1  for 
all  g  €  9,  Q  is  called  an  aitolvie  invariant  We  should  add  that  the  development  of  invariants  of  viewing 
transformations  (Euclidean,  similarity,  aASne,  and  projective  maps)  has  received  extended  attention  from  the 
computer  vision  community  in  the  last  years.  The  topic  is  very  important  in  areas  such  as  object  recognition. 
A  special  class  of  invariants  is  given  by  the  differenlittl  invariants  which  are  based  only  on  local  information, 
and  are  useful  for  the  representation  and  recognition  of  objects  under  partial  occlusions.  An  example  of 
this  kind  of  invariant  is  the  Euclidean  curvature.  In  our  work,  we  have  been  especially  interested  in  affine 
differential  invariants. 

The  third  component  of  our  work  comes  from  the  theory  of  plane  curve  evolution.  Here  the  curve  regarded 
as  the  boundary  of  a  planar  dom^,  deforms  in  time.  This  deformation  is  governed  by  a  partial,  usually 
nonlinear,  evolution  equation.  Different  evolution  equations  can  model  different  physical  phenomena,  such 
as  crystal  growth,  the  Huygens  principle,  and  curve  shortening  processes.  The  theory  has  been  well  studied 
in  areas  such  as  computational  physics,  differential  geometry,  numerical  analysis,  and  parabolic  equations. 
We  have  introduced  these  ideas  into  computer  vision  in  [49,  50,  51]. 

Motivated  by  the  importance  of  affine  transformations  in  computer  vision,  a  new  theory  of  aASne  curve 
evolution  was  recently  developed  [60,  61,  62,  63].  This  theory,  based  on  aASne  differential  geometry  and  the 
theory  of  parabolic  evolution  equations,  constitutes  the  basis  of  the  work  on  multi-scale  representations  of 
rignals  and  image  processing  in  general.  This  is  based  in  turn  on  a  new  nonlinear  evolution  diffusion  equation 
which  admits  aASne  invariant  solutions.  The  solution  of  this  evolution  equation  determines  an  affine  invariant 
scale-space  for  planar  curves.  Efficient  computer  implementation  of  this  theory  is  possible  due  to  a  recently 
developed  numerical  algorithm  for  curve  and  surface  evolution.  In  our  continuing  ARO  research,  we  will  be 
emplcqring  these  algorithms  to  a  number  of  problems  in  active  vision,  and  visual  control. 
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Analysis  and  Applications. 

16.  “On  aflSne  plane  curve  evolution”  (with  G.  Sapiro),  Journal  of  Functional  Analysis  119  (1994),  79-120. 

17.  “AflSne  invariant  scale-space”  (with  G.  Sapiro),  International  Journal  of  Computer  Vision  11  (1993), 
25-44. 

18.  “Causality  in  commutant  lifting  theory”  (with  C.  Foias),  Journal  of  Functonal  Analysis  118  (1993), 
407-441. 
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19.  "On  invariant  curve  evolution  and  image  analyeia”  (with  G.  Sapiio),  Indiana  Univ.  Journal  of  Matke- 
maiia  42  (1993),  985-1009. 

20.  "Nonlinear  IT""  optimifation:  a  causal  power  series  approach"  (with  C.  Foias  and  C.  Gu),  to  appear  in 
SIAM  J.  Control  and  OpUmvuMon. 

21.  G.  Sapiro  and  A.  Tannenbaum,  "Area  and  length  preserving  geometric  invariant  scale-sapces,”  to  appear 
in  IEEE  Trano.  on  Pattern  Analysis  and  Machine  Intelligence. 

22.  "On  a  lifting  theorem  for  the  structured  singular  value”  (with  H.  Bercovici,  C.  Foias,  and  P.  Khargonekar), 
to  ^pear  in  Journal  of  Math.  Analysis  and  AppUeatione. 

23.  "A  relative  Toepliti-Hausdorff  theorem*  (with  H.  Bercovici  and  C.  Foias),  to  appear  in  in  Operator 
Theorg:  Advances  and  Appieetione. 

24.  "Continuity  of  the  spectrum  on  closed  similarity  orbits”  (with  H.  Bercovici  and  C.  Foias),  Integral 
EguaUone  and  Operator  Theorg  18  (1994),  242-246. 

25.  "Differential  invariant  signatures  and  flows  in  computer  vision:  a  symmetry  group  approach”  (with  P. 
Olver  and  G.  Sapiro),  to  appear  as  a  book  chapter  in  Geometric  Diffueion,  edited  by  Bart  Romeny. 

26.  "Clasrification  and  uniqueness  of  invariant  geometric  flows”  (with  P.  Olver  and  G.  Sapiro),  to  appear  in 
Comptee  Rendue. 

27.  "On  optimal  control  methods  in  computer  vision  and  image  processing”  (with  B.  Kimia  and  S.  Zucker), 
to  appear  as  a  book  chapter  in  Geometric  Diffusion^  edited  by  Bart  Romeny. 

29.  "On  the  nonlinear  standard  problem”  (with  C.  Foias  and  C.  Gu),  submitted  for  publication  in  Int. 
J.  Nonlinear  and  Rohust  Control. 

30.  "Invariant  theory  and  families  of  systems,”  Mathematical  Sgttemt  Theorg:  The  Influence  of  Rudolf 
Kalman,  edited  by  A.  C.  Antoulas,  Springer- Verlag,  1991,  327-345. 

31.  "Entropy  scale-space”  (with  B.  Kimia  and  S.  Zucker),  in  Vtsval  Form,  edited  by  C.  Arcelli,  pages 
333-344,  Plenum  Press,  New  York,  1992. 

32.  "Frequency  domw  methods  for  the  ””  optimisation  of  distributed  systems,”  in  Analysis  and  Opti¬ 
misation  of  Sgeteme:  State  and  Frequencg  Domain  Approaches  for  Infinite  Dimensional  Sgstems,  edited  by 
Ruth  Curtmn,  Lecture  Notes  in  Control  and  Information  Sciences  185,  Springer- Verlag,  New  York,  1992. 

33.  "Generalised  interpolation  theory  and  its  application  to  robust  system  design,”  Digital  and  Numeric 
Technigues  and  Their  Applications  in  Control  Sgstems,  edited  by  C.  T.  Leondes,  Academic  Press,  1993, 
pages  163-217. 

34.  "From  curve  detection  to  shape  description”  (with  A.  Dobbins,  L.  Iverson,  B.  Kimia,  and  S.  Zucker), 
in  Computer  Vision:  Sgstems,  Theorg,  and  Applications,  edited  by  A.  Basu  and  X.  Li,  World  Scientiflc, 
Singapore,  1993,  pages  25-39. 

35.  "Stability  margin  optimisation  via  interpolation  and  conformal  mappings”  (with  J.  Cockbum  and  Y. 
Sidar),  submitted  for  pubUcation  in  IEEE  Tlmns.  A»t.  Control. 

Refereed  Conference  Publications 

36.  optimal  controllers  for  as  distributed  model  of  an  aircraft”  (with  D.  Enns  and  H.  Osbay),  Proceed¬ 
ings  of  SOih  IEEE  Conference  on  Decision  and  Control,  Brighton,  England,  1991. 

37.  "Computational  methods  for  the  control  of  distributed  systems,”  Proceedings  of  SOth  IEEE  Con¬ 
ference  on  Decision  and  Control,  Brighton,  England,  1991. 

38.  "On  the  optimal  two  block  compensators  for  distributed  unstable  plants”  (with  H.  Osbay),  Pro¬ 
ceedings  of  ACC,  1992. 

39.  "Conservation  laws  and  the  evolution  of  shape”  (with  B.  Kimia  and  S.  Zucker),  Human/Machine  Vision 
Workshop,  SPIE  Meetings,  Falmouth,  MA,  1991. 
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40.  "A  conaervstion-minded  approach  to  ihape”  (with  B.  Kimia  and  S.  Zucker),  Proceeiing$  of  IEEE 
Conference  on  IntelKgent  Control,  Philadelphia,  1901. 

41.  *A  constructive  solution  to  the  gw-phase  margin  problem”  (with  Juan  C!ockburn  and  Yariv  Sidar), 
Proeeeiinge  of  Slot  IEEE  Conference  on  Decision  and  Control,  Arisona,  1992. 

42.  ^Synthesis  methods  for  robust  nonlinear  control”  (with  D.  Bugajski  and  Dale  Enns),  Proceedings  of 
ACC,  1993. 

43.  "Preliminary  mu-synthesis  design  for  the  ATB-1000”  (with  D.  Bugajski  and  Dale  Enns),  Proceedings  of 
Tenth  Armg  Conference  on  AfpUed  Maihemniics  and  Computing,  West  Point,  New  York,  1992. 

44.  "Exploring  the  shape  manifold:  the  role  of  conversation  laws”  (with  B.  Kimia  and  S.  Zucker),  Proceedings 
of  the  Shape  in  Picture  NATO  Conference,  Driebergen,  Netherlands,  September,  1992. 

45.  "Causal  commutant  lifting  in  J7^(D")  and  its  applications  in  nonlinear  control”  (with  C.  Foias  and 
C.  Gu),  Proceedings  of  Workshop  on  Robust  Control,  IMA,  Minneapolis,  Minnesota,  September  1992. 

46.  "Non-linear  shape  appronmation  via  the  entropy  scale  space”  (with  B.  Kimia  and  S.  Zucker),  Proceedings 
of  the  SPIE  Geometric  Methods  of  Computer  Vision  Conference,  San  Diego,  1993. 

47.  "Formulating  invariant  heat-type  curve  flows”  (with  G.  Sapiro),  Proceedings  of  the  SPIE  Geometric 
Methods  of  Computer  Vision  Conference,  San  Diego,  1993. 

48.  "A  lifting  technique  for  the  robust  stability  analysis  of  systems  with  structured  time-varying  pertur¬ 
bations”  (vrith  H.  Bercovici,  C.  Foias,  and  P.  Khargonekar),  Proceedings  of  the  Conference  on  Information 
Sciences  and  Sgstems,  Johns  Hopkins  University,  March  1993. 

49.  "Image  smoothing  baaed  on  an  afline  invariant  flow”  (with  G.  Sapiro),  Proceedings  of  the  Conference  on 
Information  Sciences  and  Sgstems,  Johns  Hopkins  University,  March  1993. 

50.  "On  structured  tangential  interpolation  in  robust  control”  (with  H.  Bercovici,  J.  Cockbum,  and  C. 
Foias),  Proceedings  of  Stnd  IEEE  Conference  on  Decision  and  Control,  December  1993. 

51.  "Pseudorational  functions  and  H*”  theory”  (with  Y.  Yamamoto),  to  appear  in  Proceedings  of  ACC, 
1994. 

52.  "Skew  Toeplits  theory  and  pseudorational  transfer  functions”  (with  Y.  Yamamoto),  submitted  to  1994 
IEEE  CDC. 

53.  "Experiments  on  geometric  image  enhancement”  (with  M.  Kaveh,  G.  Sapiro,  Y.  L.  You),  submitted  to 
First  IEEE  International  Conference  on  Image  Processing,  Austin,  1994 

Books 

1.  J.  Doyle,  B.  Fkancis,  and  A.  Thnnenbaum,  Feedback  Control  Theory,  MacMillan,  New  York,  1991. 

2.  C.  Foias,  H.  Osbay,  and  A.  'Ihnnenbaum,  Control  of  Distributed  Parameter  Systems,  in  preparation, 

1994.  First  draft  is  avulable. 

8  Students  of  Allen  Tannenbaum  Supported  by  DAAL-03-91-G- 
0019 

1.  Juan  Cockbum  (Ph.  D.) 

2.  Arun  Kumar  (Ph.  D.) 
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